Levy statistics and anomalous transport in quantum dot arrays 
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A novel model of transport is proposed to explain power law current transients and memory phe- 
nomena observed in partially ordered arrays of semiconducting nanocrystals. The model describes 
electron transport by a stationary Levy process of transmission events and thereby requires no time 
dependence of system properties. The waiting time distribution with a characteristic long tail gives 
rise to a nonstationary response in the presence of a voltage pulse. We report on noise measure- 
ments that agree well with the predicted non-Poissonian fluctuations in current, and discuss possible 
mechanisms leading to this behavior. 

PACS numbers: 73.50.Fq,73.61.Ga,73.63.Kv,73.50.Td 



I. INTRODUCTION 

Arrays of semiconductor nanocrystals 1 (quantum dot ar- 
rays, or QDAs) are of great interest for both the fun- 
damental solid state physics and applications. Self- 
assembled QDAs is one of the simplest examples of 
macroscopic complex systems built from the "artificial 
atoms" with pre-designed properties at the nanoscale. 
From the basic research perspective, these arrays are 
compelling since one can control the Hamiltonian by de- 
sign. In particular, they open new possibilities to create 
systems with desirable unconventional transport proper- 
ties. Charge and spin transport in QDAs could lead to 
applications in spintronics and quantum computation. 2 

Despite the progress in synthesis and fabrication of 
nanocrystal arrays, the nature of electronic transport in 
them is still poorly understood. Like other problems in- 
volving long-range Coulomb interactions in disordered 
systems, this problem appears to be challenging. Pro- 
posed theoretical models include mapping onto the prob- 
lem of interface dynamics, 3 onto a frustrated antiferro- 
magnetic spin model with long range interactions, 4 as 
well as recent generalizations 5 of the variable range hop- 
ping scenario. 6 

In the present work we attempt to explain the recently 
observed 7 ' 8 transient power-law decay of current 



/(f) = i r 



< a < 1 



(1) 



as a response to a step in large bias voltage applied across 
the array. The exponent a depends on temperature, dot 
size, capping layer, bias voltage and gate oxide thick- 
ness in a systematic way. 7 ' 8 It is interesting that the 
observed a is less than one in all samples. The con- 
dition a < 1 ensures that (1) is a true current from 
source to drain, rather than a displacement current, since 
the net charge corresponding to (1) diverges with time, 
Q = J I(t)dt oo. At the same time, remarkably, the 
transient transport also possesses memory. Namely, if 
the bias is turned off for t\ < t < t^, then the current 
measured as a function of the time t = t — ti after reap- 



plying the voltage follows the dependence (1), albeit with 
a reduced amplitude Io < Io- I{t) = I (t — t2)~ a - The 
amplitude is restored to its initial value, Io — » Io, by in- 
creasing the off interval £2 — t\ , by annealing at elevated 
temperature, or by applying a reverse bias or band gap 
light between t\ and ti. 

The behavior (1) is observed in partially ordered 
multi-layered arrays of II- VI semiconductor nanocrys- 
tals. Each nanocrystal is capped with ~lnm coating, 
so that electrons must tunnel to move between neighbor- 
ing sites. Although the transient response (1) has been 
recorded by a number of groups, 7-11 its origin remains 
a mystery. In many systems the transient (1) comprises 
the dominant contribution to transport, while the ohmic 
conductivity is quite small. 12 Understanding the nature 
of the time-dependent response (1) may shed light on the 
dynamics of carriers in such systems. 

It has been suggested earlier that the observed time- 
dependent current could be a result of time dependence 
of the state of the system. The latter could arise either 
because of charge flow jamming, due to trapping of elec- 
trons blocking further charge injection from the contact, 9 
or because of the Coulomb glass behavior of the electrons 
distributed over QDAs. 7 However, it seems that such a 
scenario would require an unlikely fine-tuning as a func- 
tion of time. Namely, the system's properties would have 
to adjust in a coherent fashion over many hours to yield 
well-reproducible power laws in current, observed over at 
least five orders of magnitude in time in a broad variety 
of samples. 

The purpose of this work is to suggest an alternative 
point of view on transport in QDAs, which does not re- 
quire time-varying system's properties. We propose a 
model, based on the Levy statistics of waiting times be- 
tween charge transmission events, in which the system 
remains stationary in statistical sense, but nonetheless 
exhibits a transient response. The model is corroborated 
by the measurements of the spectrum of time-dependent 
current fluctuations in CdSe QDAs, and a good agree- 
ment is demonstrated. 
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This paper is organized as follows. We begin with in- 
troducing a phenomenological model of transport which 
yields the response (1). We consider manifestations of 
this transport mechanism in the noise spectrum, and 
report the results of noise measurements. Finally, we 
briefly discuss possible microscopic mechanisms that 
could be consistent with our transport model. 



II. MODEL OF TRANSPORT 

The main idea of our approach is that current (1) can 
arise as a result of a stationary stochastic process. Our 
model involves N 3> 1 identical independent conducting 
channels arranged in parallel. (This accounts for the typ- 
ical sample's aspect ratio ~ 10 3 /Ltm: 1 /xm.) Each chan- 
nel is almost always closed, and opens up at random for 
a short interval r to conduct a current pulse that cor- 
responds to a unit transmitted charge, as schematically 
shown in the lower inset of Fig. 1. We further assume 
that the intervals between subsequent transmissions are 
uncorrelated, making the process completely character- 
ized by the waiting time distribution (WTD) of intervals 
between successive pulses. 

In particular, we will be interested in WTD p{r) with 
a broad tail at long times. In order to model the power 
law decay of the current transient, here we consider a 
special form of WTD with a long tail of the Levy type: 

p(r » r ) ~ , < M < 1 , (2) 

with To a miscroscopic time scale. Note that all moments 
of p(r) diverge. The behavior of the WTD at short times, 
p(r ~ to), is not of interest, since it does not affect the 
long time dynamics. 

As shown in Appendix A, WTD of the form (2) in- 
deed yields the power law decay (1) for the mean value 
of the current in a single channel with a = 1 — fi and 
Iq = fisimrfi/iTa. Qualitatively, the decrease in current 
with time can be understood as follows. The mean value 
of the waiting time for the process with WTD (2) is in- 
finite. Thus, for a stochastic process governed by the 
WTD (2) which started infinitely early in the past, the 
observed value of the current would be zero. Turning the 
bias on at t = sets the clock for the process. In this 
case, for the measurement interval t, only the waiting 
times t < t can occur, as illustrated by the simulation 
shown in Fig. 1 (note the double log scale). Observing 
the current over a larger time period effectively increases 
the chances for a channel to be closed for a longer time in- 
terval, which results in the decay in the average current, 
the latter approaching zero at t — > oo. We note that in 
this transport model the system's parameters character- 
izing the distribution (2) are time independent. Hence 
the process is stationary, i.e. each charge transmission 
event occurs after a delay time t described by the distri- 
bution p(t) independent of the total time t passed after 
the beginning of the measurement. 




time f, a.u. 



FIG. 1: Time dependence of the net transmitted charge Q(t) 
in a single channel (red line) and charge (Q{t)) averaged over 
N — 100 channels (blue line) simulated using WTD of the 
form (2) with — 0.5 (double log scale). Dashed line is 
a power law Q oc t M . Upper inset: The same plot in the 
linear scale. The large charge noise in a single channel is 
due to the lack of self-averaging for a wide-tail WTD. Lower 
inset: Current in a single channel with a wide distribution 
of waiting times (schematic) . Short packets of current pulses 
are separated by very long waiting times r > ro. 



Continuous time random walks with the Levy WTDs 
often arise in the systems characterized by wide distri- 
butions of time scales. 13 In the semiconductor physics 
the Levy processes have been extensively studied in 
the context of the dispersive transport, e.g. in amor- 
phous semiconductors. 14 A simple example is a system 
of electrons moving between charge traps. Its dynam- 
ics depends on energy e via an activation exponential, 
t = Toe^ £ , where (3 is the inverse temperature. For the 
distribution of the energies e described by the density of 
states of exponential form, v(e) = vqc^^ , one obtains 
p(t) of the power law form (2) with exponent /i = b/ (3 
and a = vq Tq //3. 

The probability distribution (2) leads to an unusual be- 
havior which is the subject of the theory of Levy flights. 15 
The main characteristic of the Levy statistics 13,15 is the 
violation of the central limit theorem. To illustrate this 
unconventional behavior, let us recall what happens in 
a Poissonian channel characterized by the finite mean 
waiting time f . The mean value of the transmitted 
charge Q grows linearly with time, iff) = t/f, corre- 
sponding to a constant current. The variance of the 
charge is proportional to the mean, ((Q 2 )) = \ (Q), 
in other words the relative charge fluctuation decreases, 

{{Q 2 )) 1 ^ 2 I {Q) t~ x l 2 , in accord with the central limit 
theorem. In contrast, in the case of the distribution 
(2), the mean transmitted charge increases sublinearly 
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as (Q) ~ t^, whereas its variance is proportional to the 
square of the mean, ((Q 2 )) oc (Q) 2 (see Appendix A). 
Since relative charge fluctuation does not decrease with 
time, transport in a single channel is dominated by large 
fluctuations of waiting times (see Fig. 1). 

Although in our model any given channel lacks self- 
averaging, the charge summed over N 3> 1 independent 
parallel channels averages to a smooth power law, Fig. 1, 
with fluctuations reduced by a factor of N~ x / 2 . For the 
typical sample geometry used in our experiments, con- 
sisting of ca. 50 layers, each layer of 1.6 • 10 5 dots wide 
and 200 dots across, one expects large effective N and 
small current fluctuations, as in Rcfs. 7,8. 

Also, as a check of robustness of this scenario with 
respect to spatially varying system properties, we con- 
sidered parallel non-identical channels, characterized by 
non-equal values of \i. We found that the average cur- 
rent obtained from such a model is approximately de- 
scribed by a power law of the form (1). We performed 
a simulation with TV = 100 channels, with the expo- 
nents of different channels drawn from a flat distribution, 
0.45 < \i < 0.55. In this case, the transmitted charge 
time dependence was found to be numerically very close 
to that given in Fig. 1 with \i = 0.5. 

III. MEMORY EFFECTS 

Memory effects originate in our transport model in the 
manner analogous to aging in the Levy systems. 16 Be- 
cause of large typical waiting times r 3> t , any given 
channel is most likely found in a non-conducting state 
when the voltage is turned off at t = t\. In addition, 
we assume that, due to gradual time variation of channel 
parameters, taken to be very slow in this discussion, the 
channel state is likely to remain unchanged by the time 
the voltage is turned back on at t = t 2 >t\. In this case 
the channel conducts current as if the voltage has been on 
all the time. However, due to the aforementioned time 
variation of system parameters, there is a chance that 
the channel changes its state (resets) while the voltage 
is turned off during t\ < t < t 2 . This reset probability 
Wi2 = w(t 2 — t\), as a function of the off time t 2 — t\, 
is growing monotonically: w(t ) w 0, w(oo) = 1. As a 
simple model of this behavior, one can consider a Poisson 
process, 

w(t) = 1 - e" 7 * 

with the rate parameter 7 characterizing the reset prob- 
ability. 

The current at t = t 2 as a function of the shifted time 
i = t — t 2 , obtained by averaging over N 3> 1 channels, 
is given by 

7(f) = (l-w 12 )I (t + t 2 )- a +w 12 I r a . (3) 

Here we assume the reset of different channels to be in- 
dependent and uncorrelated. The function I(t) has a 



singular part at t « (the second term of Eq.(3)) with 
the amplitude 7o = wi2 7o reduced compared to 7o by 
the reset probability w\ 2 < 1. Thus at t <C t 2 — t\ the 
first (regular) term in Eq. (3) is negligible compared to 
the second term. The current (3) is dominated by the 
latter, resulting in suppression of the measured transient 
current part which is singular at f w t 2 . 

We note that the described reset process, while lead- 
ing to suppression of the singular part of the current, is 
accompanied by an overall enhancement of the total cur- 
rent (3), as compared to the current (1) at time t in the 
absense of resetting. This prediction indeed agrees with 
our observations. We have verified that the reset proba- 
bility w\ 2 = 7 /7 is indeed a monotonic function of the 
time interval when the voltage is turned off. For waiting 
times from 10 s to 10 4 s in between 100 s long transients, 
we measure 0.65 < w\ 2 < 0.85; wi 2 —> 1 when applying 
a reverse bias, exposing the dots to the band gap light or 
waiting for longer times. 

IV. NOISE FREQUENCY SPECTRUM 

The model described above, which is consistent with pre- 
viously reported transport measurements, can be inde- 
pendently verified with the help of noise measurements. 
Here we consider the statistics of current fluctuations and 
formulate a prediction of the model based on the Levy 
process (2). 

The unconventional fluctuations exhibited by the Levy 
process, discussed in Section II, manifest themselves in 
noise as follows. Consider the time-dependent current in 
a single channel, described in Section II, recorded during 
a long time interval T » r : 

I(t)= S (t~tn), 0<t<T. (4) 

71=1,2,... 

The intervals r„ — t n — f n _i, n = 1,2,..., to = 0, are 
independent random variables distributed according to 
the WTD of the form (2). The fluctuations of current 
are defined in terms of the Fourier harmonics 

I u = [ e wt I{t)dt . (5) 
Jo 

Here we consider the noise power spectrum 

S(u) = «7_ W 7 W )) = (7^7,) - (J_ w ) (/„) . (6) 

In Appendix A we show that the distribution (2) leads 
to the non-Poissonian spectrum 

Here the low frequency part of (7) with ojT <C 1 corre- 
sponds to the fluctuation ((Q 2 )) of the net transmitted 

charge. Due to the relation ((Q 2 )) oc (Q{T)) 2 [Eq. (A14)] 
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between the first two moments of the Levy process, 
which violates the central limit theorem, the quantity 
S(ujT <C 1) is proportional to the square of the mean 
transmitted charge Q(T) oc T M . 

Experimentally, however, it is more convenient to deal 
with S(ui) at finite frequency ujT 3> 1. Eq. (7) predicts a 
characteristic power law spectrum for this quantity. We 
note that for ojT ;§> 1, the Levy process (2) yields iden- 
tical power laws for the noise spectrum (7) and for the 



average current, (I J) 



l . Moreover, the relationship 



(8) 



is robust with respect to averaging over TV independent 
channels, since for such averaging the central limit the- 
orem holds. An experimental test of the proportionality 
relationship (8) between the frequency spectra of current 
and noise will be discussed below. 



V. NOISE MEASUREMENTS 

Here we briefly describe the experiments performed to 
obtain the data on noise frequency spectrum in QDAs. 
The QDAs were produced as described in Ref.7 by self- 
assembly of nearly identical CdSe nanocrystals, 3 nm in 
diameter, capped with trioctylphosphinc oxide, an or- 
ganic molecule about 1 nm long. A film of about 200 
nm thick of the nanocrystals was deposited on oxidized, 
degenerately doped Si wafers with oxide thickness ks 
200 nm. The experimental setup was similar to that uti- 
lized in Ref. 7. Gold electrodes, fabricated on the sur- 
face before deposition of the QDA, consist of bars 800/mi 
long with separation of 2/zm. The sample was annealed at 
300 C in vacuum inside the cryostat prior to the electrical 
measurements. Annealing reduces the distance between 
the nanocrystals and enhances electron tunneling. 7 

To measure the noise, we have recorded 200 current 
transients each t = 100 s long. Measurements have been 
made on a single sample continuously stored in vacuum, 
inside of a vacuum cryostat in the dark at 77 K. Each 
current transient was recorded for 100 s with a negative 
bias of —90 V. These periods of negative bias were sepa- 
rated from each other by a sequence of zero bias for 10 s, 
reverse pulse of +90 V for 100 s, and zero bias for 10 s, 
to eliminate the memory effects. 7 We checked that cur- 
rent fluctuations for a substrate without the QDA were 
several orders of magnitude smaller than with the QDA. 

At the beginning of our measurement the current tran- 
sients were changing from one to the next because of the 
memory effect described above. Since an error in the av- 
erage current can yield an error of order {I u ) 2 oc ui~ 2fl 
which may affect the measured noise spectrum power 
law, we discarded the first 150 transients. The noise 
(Fig. 2) was then deduced from the remaining 50 tran- 
sients. To further compensate for residual memory ef- 
fects, each transient was multiplied by a factor rs 1 to 
have the same net integrated charge. This eliminated 
the zero frequency contribution to the noise. 




Frequency f=G>/2iz, Hz 



FIG. 2: Measured mean current and noise spectrum, with 
averaging performed over 50 current transients on the same 
sample. The current Fourier harmonic mean [Eq. (5)] 

and variance ((1-^1 J)) [Eq. (6)] are shown. Both the current 
and the noise are described by power law with the same 
/i w 0.72. The 1/f dependence is shown for comparison (see 
text). 



Figure 2 shows the measured noise spectrum and the 
average current measured simultaneously. Both quanti- 
ties have a power law behavior with nearly identical expo- 
nent values /i ~ 0.72 for u>t 3> 1, t = 100 s. With relative 
~ 3 — 10% deviations of the noise from the law, the 
measured noise spectrum is clearly distinct from the 1 // 
noise typically found at low frequencies. For compari- 
son, in Fig. 2 we draw the 1/f dependence, offset so that 
it coincides with the noise data at the lowest frequency. 
The discrepancy with the measured noise at the highest 
frequency by more than a factor of two indicates that the 
observations are not explained by the 1/f noise model. 

One may question whether the observed colored noise, 
instead of being a consequence of the Levy process, could 
result from interplay of the intrinsic 1// noise and the 
time-dependent current decaying according to (1). In 
this case the fluctuations would be proportional to the 
current itself: 

I(t) = {I(t))(l + s(t)), (8{tW)) = T{t-H), (9) 

where r w oc uj^ 11 , 77 ~ 1 for the 1/f noise. This would 
yield the current fluctuation spectrum of the form 

«/-„/*» = / \{i*>)\ 2 r„,-^. (10) 

When /i > 1/2, the integral is dominated by the (a/) -2 ^ 
singularity, which effectively sets u>' = 0, giving rise to 
the 1/oj behavior. Physically this happens because the 
current (1) for a < 1/2 decays slowly enough so that all 
the harmonics of the 1/f noise have time to fully play 
out. 
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Conversely, in a system with /i < 1/2, the noise oc io^ 2 ^ 
may be indistinguishable from the errors in determin- 
ing average current described above. For the consistency 
check of our model it is important that for our sample 
the observed current transient power law exponent fulfills 
[i > 0.5. The observation of the uj~^ noise indicates that 
transport is not dominated by the intrinsic 1// noise. 
Instead, we conclude that the noise measurement agrees 
with the proposed transport model based on Levy statis- 
tics of transmission events. 



VI. DISCUSSION 

Based on the noise measurements we can estimate the 
effective number TV of conduction channels introduced 
in Section II. Since both the measured current (J w ) and 
noise S(lo) are proportional to N, the noise-to-current ra- 
tio r{uo) = S(oj)/ is a characteristic of a single chan- 
nel. The model calculation for the latter [Appendix A] 
shows that at large frequency ut 3> 1 the noise-to-current 
ratio is frequency independent and proportional to the 
net charge Q(t) transmitted through the channel during 
the time t of measurement, r(wt 3> 1) ~ Q(t) [Eq. (A15)]. 
From Fig. 2 we find that the measured r(u>) is indeed 
frequency independent for > 1, and is of the order 
r w 10 -10 A/Hz. [Averaging over 50 transients does not 
affect r(w).] The effective number of channels is then 
estimated as the ratio of the measured net transmitted 
charge to that in a single channel, N ~ (I w= o) /r ~ 10 4 . 
This large number of independent channels is consistent 
with the sample geometry (aspect ratio ~ 10 3 and ca. 50 
layers of dots). 

How can the long waiting times with a distribution 
of Levy form (2) arise microscopically? While presently 
there is no fully satisfactory answer to this question, 
one can make several observations. First, to rationalize 
a wide distribution of time scales such as (2), we sup- 
pose that the charge hops between neighboring dots are 
strongly constrained. The simplest constraint to imagine 
is the lack of energy relaxation (possibly due to a small 
number of available phonon states) that arises if the on- 
site energies of electrons on different dots are widely dis- 
tributed. This is consistent with the absence of ohmic 
contribution to the conductivity in our QDA. The en- 
ergy relaxation constraint then allows charge hops only 
between the aligned energy levels of the dots. 

Next, the WTD (2) with a long tail can be explained if 
the energy levels strongly fluctuate in time, with [i. = \ 
corresponding to the Gaussian diffusion in energy. One 
can think of at least two reasons for the level fluctua- 
tions. First, the voltage bias energy ~ 0.1 eV dissipated 
per hop may provide the necessary energy reservoir. Sec- 
ond, current induced fluctuations in the electrostatic en- 
vironment in the absence of screening may result in a 
random time-dependent chemical potential for each dot. 
In particular, misalignment of the energy levels can arise 
due to the Coulomb field of an electron trapped in the 



vicinity, e.g. in the coatings. The "conduction channel" 
then opens up when the trapped electron escapes. Due 
to large applied bias, filling of the traps can happen much 
faster than escaping from them. With escape times expo- 
nentially dependent on trap parameters, the distribution 
(2) follows naturally. 14 

We note that this picture differs somewhat from the 
canonical dispersive transport mechanism, 14 in which a 
constant supply of carriers makes the current grow with 
time. 17 The growth of current occurs due to the increas- 
ing number of trapped electrons levelling the potential 
landscape, thereby enhancing conductivity. Contrarily, 
in the proposed picture the presence of traps regulates 
the dynamics of conducting channels. 

We also note that the Levy statistics was recently 
observed in fluorescence intermittency of individual 
nanocrystals. 18 Possibly, a better understanding of the 
microscopic mechanism of the anomalous transport can 
be achieved by establishing a connection between the 
statistics of fluorescence and of charge transmission in 
the same sample. This could discriminate between trans- 
port due to the properties of electron states in a single 
nanocrystal, and the collective transport phenomena. 



VII. CONCLUSIONS 

This article presents a novel mechanism for a non-ohmic 
conductivity in a disordered system. In particular, we 
show that a non-stationary current response can arise 
in a stationary system with the Levy statistics of wait- 
ing times. The model agrees well with the current and 
noise measurements in arrays of coated semiconducting 
nanocrystals. The non-Poissonian character of the Levy 
process manifests itself in the non-ohmic character of 
transport observed as the current transients, in mem- 
ory effects, and in the colored noise. Our results suggest 
that one needs to be careful in interpreting conductivity 
in such systems 10 using simple ohmic models implying 
Poissonian statistics of transmission. We also demon- 
strate that the Levy model can help to investigate the 
system even without precise knowledge of microscopic 
transport mechanism, by linking the power law observed 
in the noise with that of current transient. 
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APPENDIX A: CURRENT AND NOISE IN A 
SINGLE CHANNEL 



The Fourier harmonic variance, obtained from Eq. (A6), 
is of the form 



Consider the current in a single channel 

oo 



(Al) 



71=1 



Here, instead of switching the current on and off at t = 
0, T as in (4), we introduced a soft cutoff X = T -1 . This 
cutoff helps to simplify calculations without qualitatively 
affecting the results. The Fourier harmonic of (Al) is 



oc 

Iu = j e-^I(t)dt = J2' 



(A2) 



where z = A — itu and t n = Y^i=i Ti • Since waiting times 
n are independent random variables distributed accord- 
ing to p(r), average current is given by the geometric 
series 



Pz 



with p z the characteristic function 



Pz = (e- ZT ) = J e- ZT p(r)dr . 



(A3) 



(A4) 



The correlator (J_ W J W ) = E^„'=i ( e_Zt "' _z '") can be 
evaluated as 



rn—1 



(A5) 

The last formula is obtained by splitting the summation 
into parts with n = n' and n > n' with m = n — n' . The 
variance = (I- U I U ) - {I-J) is given by 



((/ i )) - P2X-PzP-z 

{{ UU)) (1-P2A)(1-P,)(1-P 2 )' 



(A6) 



where p 2 \ = p z =2\- 

The expressions for current average (A3) and variance 
(A6) are valid for any waiting time distribution. Con- 
sider now the WTD of the form (2). In this case the 
characteristic function is 



Pz = 1 - -7- , A„ 



aT(l - n) 



(A7) 



where \z\t <C 1, corresponding to the long time tail. 
Eq. (A3) then yields 

(/„> = (A — iw)-", (A8) 

resulting in the average current of the form (1): 

(I(t))=I r a , a = l-n, 1 = ^^. (A9) 

na 



((/_„/„» = Al 



- {2\y - A~ 1 \z\ 2 ^ 

{2\y \z\ 2 ^ 



. (A10) 



We are interested in the noise spectrum on the time scale 
much greater than the pulse width t$: uj, A <C Tq 1 . Keep- 
ing the leading terms in Eq. (A10), we have 



ZM + ZM_(2A)M 



(All) 



The limits of Eq. (All) arc [T = A" 1 ]: 

nr r u_/(2 1 -''-l)4 T2 ' i . wT « 1 ' 

(A12) 

The result (A9) corresponds to the net transmitted 
charge 



(Q) = A^ . 



(A13) 



For the exponent p < 1, all the moments of (2), including 
the mean and the variance, diverge, and thus the central 
limit theorem does not hold. Instead of (Q) oc T expected 
for a stationary random process, here we have a power 
law. Moreover, the distribution of Q(T) is extremely 
broad, with dispersion proportional to the net charge: 



varQ- ((Q 2 )-(QY 



1/2 



(2- 



i) 1/2 (Q) ■ 



(A14) 

Hence the ratio varQ/ (Q) does not decrease with time 
T, violating the central limit theorem. 

The large frequency asymptotic behavior of the cur- 
rent and noise is given by the same characteristic power 
law. According to Eq. (A12), the noise-to-current ratio is 
controlled by the net transmitted charge (A13) through 
the channel: 



= fc^=cQ ( T), C = 2^cosf . (A15) 



It is instructive to compare our results for the WTD 
(2) with those derived for the Poissonian statistics. In 
the Poissonian case, p(r) — f~ 1 e~ T / T , p z = (1 + zf) _1 . 
Eq. (A3) yields the average current {I J) = [(A — 2w)f] _1 , 
corresponding to (I(t)) = 1/f for t < T = A -1 . At 
the same time, Eq. (A6) yields the white-noise spectrum 
((I-ulu)) = T/2f = 2 (Ilj=o), in agreement with the 
central limit theorem. 
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